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CHART FOR CRITICAL COMPRESSIVE STRESS OF 



FLAT RECTANGULAR PLATES 



By H. N. Hill 
SUMMARY 



A chart is presented for the coefficient Z in the 
formula for the critical compressive stress for flat rec- 
tangular plates uniformly compressed in one direction. 
The chart applies to various combinations of fixed, simply 
supported, and free edges. 



Chart for K in Formula for Critical Stress 



The theoretical stress at which a rectangular flat 
plate will "buckle elastically, when uniformly compressed 
in one direction, can "be expressed (see reference 1, p. 
331) 



cr cr = E 



il-u 3 ) W 



where <j cv critical stress (at which buckling occurs) , 

pounds per square inch 

E modulus of elasticity of the material of the 
plate, pounds per square inch 

H Poisson 1 s ratio of the material of the plate 

t thickness of the plate, inch 

b width of the plate (normal to the direction 
of compression) , inch 

E a coefficient depending on the ratio of length 
to width of the plate (L/b) , the conditions 
of restraint at the edges of the plate, and 
in some cases on the Poisson's ratio of the^ 
material of the plate 
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In designing structures involving flat plates, a 
chart from which K values for the above equation could 
be obtained, for various conditions of edge restraint, 
should prove a convenience. ' It would be impossible to 
include all the possible variations or combinations of 
edge conditions, but by defining three definite edge con- 
ditions, and considering the various combinations of these 
conditions that might occur, certain limiting cases are 
obtained. Other cases involving intermediate conditions 
of edge restraint will be bounded by several of those lim- 
iting cases. The definite ■ edge conditions considered are: 

1. Free edge - an edge about which the plate is free 

to rotate and at which the plate is free to de- 
flect. 

2. Supported edge -an edge about which the plate is 

free to rotate but at which thero can be no do- 
flection. 

3. Fixed edge - an' edge at which there can bo no de- 

flection and about which thoro can bo no rota- 
tion of the plato. 

The buckling stress for a rectangular flat plate sub- 
jected to uniform compression in one direction, with the 
unloaded edges free, can be determined from the column 
formulas. For plates, "E " in the column formulas should be 
replaced by the 'plate modulus" [E/ ( 1-u- 3 ) ] (reference 2, p. 475). 

For plates with other edge conditions E values to be 
usod in the equation for critical stress can bo obtainod 
from the curves plotted in figure 1_. Curvos are givon for 
five combinations of conditions of restraint at the unloaded 
edges. for each of these combinations two conditions of 
the loaded edges have been considered. Theoretically, 
curvos for all casos but cases 1 and la should bo composed 
of a scries of sca.llops, each scallop representing buckling 
into a certain number of half waves. For practical reasons, 
however, these scallops have been eliminated. The curves 
shown in figure 1 arc tho envelope curves drawn tangent to 
the scallops. 

The curves of figure 1 represent various degrees of 
approximation to the theoretically correct values, although 
in every case the curves are sufficiently accurate for de-_. 
sign purposes. The sources from which K values for tho 
different cases were obtained are given in table I, together 
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with pertinent remarks concerning the plotting of the 
curves. For those cases for which K values were obtained 
directly "by the writer, the major elements of the method 
of solution employed are given in table I and the solutions 
are included in an appendix to this report. 

For the cases in which one unloaded edge of the plate 
is free (1, la, 2, and 2a), the K values are based on a 
Poisson's ratio value of 1/3, which is the accepted value 
for aluminum alloys. In some instances, a difference be- 
tween Poisaon's ratio values of l/4 and l/3 is responsible 
.for a difference of as much as 15 percent in the K value. 
Por cases other than those involving one free edge, the "" 
curves of figure 1 are applicable to ary elastic material, 
the K values being independent of Poisson's ratio. 

Che curves of figure 1 demonstrate the importance of 
the condition of restraint of the loaded edges in determin- 
ing the buckling stress of rectangular flat plates having 
L/b ratios lossr than- 3. 



Aluminum Company of America, 
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TABLE I 



SOURCE OP K VALUES- PLOTTED 11?" FIGURE 1 



Case 



Source and Remarks 



1 Solution from Timoshenko »'s : "Theory of Elastic 

Stability , 11 pp. 339-340, and the values calculated 
by the writer for a P'ois-son's ratio value of l/3. 

la Approximate solution by the writer, using the energy 
method and the deflection function 

w = A(l-cos 2n-x/L)(y + By?), 'in which B = f(b,n) 
is determined from boundary • condit ions . 

2 Solution from Timoshenko's "Theory of Elastic 

Stability," p. 341, and values calculated by the 
writer for a Poisson's ratio value of l/3. 



2a Approximate solution by the writer, using the energy 
method and deflection functions 

w = A(l-cos 2ttx/L) (l-cos ny/b) for buckling in one 
half wave, and w = A(sin 2Jx/L + Bx 3 /L 3 )( l-cos ny/B) 
for buckling in any even number of half waves. Coof- 



c i ent s n , 
condit ions 



and J are determined from boundary 



Solution from Timoshenko's "Theory of 
Stability," pp. 329-332. 



Elastic 



3a Solution from Timoshenko's "Theory of Elastic 
Stability," pp. 363-364. 

4 Solution by the writer following the method employed 
by Timoshenko in "Theory of Elastic Stability," 
pp. 337-342. 

4a Tho relation if this curve to that for case 4 is esti- 
mated from tho relatione between the curves for 
casos 3 and 3a and cases 5 and 5a. 



5a 



Solution from Timoshenko's "Theory of Elastic 
Stability," pp. 344-345. 

Solution from "Buckling of Compressed Rectangular 
Platos with Built-in Edges," by J. L. Maulbetech, 
Journal of Applied Mechanics, June 1937. 



HA OA Technical Soto No. 773 



5 



APPENDIX 



Consider a flat plato of length L and width. b, 
subjected to a uniform end compression (in the L direc- 
tion) of N s per unit width. Consider the axes disposod 
as indicated in figure 2. If the deflection of the plate 
out of its original plane is denoted "by w, the differen- 
tial equation for the buckled plate may he expressed (ref- 
erence 1, p. 337) 

3 4 w , 0 9 4 w , 5 4 w ^x 5 S w /, \ 

al* + 3 ax 2 a y a + a y 4 ~ ~ d aT* w 

E t 3 

where D = — rigidity of the plato per unit width 

l2(l-|i S ) 

The various boundary conditions for the edges of the 
plato parallel to the x axis may be expressed as follows 
(reference 1, pp. 293-300): 

Por '£ ixed" edges , 

w = 0 and ~ = 0 

3Por "supported" edges, 

w = 0 andf^f+^!^!=0 
SyS ^ 8x3 

Por "free" edges, 

|f| + „ |f| = o and |4 + (3-n) = 0 

By 3 Bx2 By 3 3y3x s 

In cases in which load is applied to " "supported" edges 
the problem of the buckling of a flat plate under edge com- 
pression can be solved by integration of the differential 
equation, obtaining the constants of integration from tho 
boundary conditions. Such a method has" been employed by 
Timoshenko for obtaining the solution for cases 1, 2, and 
5, as well as for certain cases involving intermediate con- 
ditions of edge restraint (reference 1, pp. 337-350). In 
obtaining a solution for case 4, the writer has followed 
the integration method of Timoshenko, determining integra- 
tion constants consistent with the boundary conditions for 
this case. 



6 



STAC A Technical Note "No. 773 



For cases in which the loaded edges are 'fixed," the 
above-mentioned method of integrating the differential 
equatdon cannot he used, and it is convenient to resort 
to one of the approximate methods for obtaining a solution. 
The writer has obtained approximate solutions for casos 
la and 2a by application of the "energy" method, This 
method also has beon used extensively by Timoshenko (ref- 
erence 1, p. 81). The "energy equation," obtained by 
equating expressions for the work of external forces and 
the energy of bonding nay bo written: 



dx dy 



/7(H) 

■Jo Jo 



(2) 



The critical value for IT X ' is obtained by solving this 
equation, using an expression for the deflection w, which 
is some function of x and y and satisfies the boundary 
conditions of the case being considered, Excopt in the very 
improbable instance when the assumed expression for the de- 
flection exactly defines the configuration of the buckled 
plate, the critical load values obtained by the energy meth- 
od are always higher than the true values (reference 1, p» 
81). The degree of approximation of tho energy method can 
sometimes bo improved by including a. parameter in the -ex- 
pression for the deflected surface, and evaluating this pa- 
rameter so as to minimize the value determined for tho 
critical load. 



Case la 

Loaded edges (x = 0 and x = L) fixed; one unloaded 
edge (y = 0) supported and the other unloaded edge (y = b) 
free. 



The 


boundary 


c ondi t ions 


for thi.s 


case 


are : 




[1] 


w = 0 


when 


x = 0 


and 


when 


x = L 


C2] 


pL = 0 
o X 


when 


x = 0 


and 


when 


x = L 


[33 


w = 0 


when 


y = 0 
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To solve this problem "by the energy method, assume 
an expression for the "buckled surface of the form, 

w = f(x) 3P(y) 

Prom boundary condition [5] 

»(x) I 11 (b) = - \it H (x) JF(b) 

(The primes indicate the order of tho partial derivatives.) 

This relation may "be written 

f (x) m , ^(0) 
f«(x) 3"'(b) 

Similarly, from "boundary condition 6 

f(x) I"'»(b) = -(2-M.) f"(x) 2"(b) 

or 

1UL. = -(2.^)11^) 
f»(x) P»'(b) 

Combining these two equations 

liS-Cb) S""(b) = {2-v) F'(b) P»(b) (s) 

Boundary conditions 3 and 4 are satisfied by an expres- 
sion of the form 

w = f(x) (y+By 11 ) 
i.e. JCy) = (y+By n ) (4) 
where n 1 
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Substituting equation (4) and its deri-vatives into equation 
(3) and solving for B 



H 



where 



H = 



2-jJin+[i 
2 (|j,n-|j.-n) 



(5) 
(6) 



The "boundary conditions [ l] and [2], defining the conditions 
at the loaded edges, are satisfied by 



f(ac) = A^l-cos ^jp) 



The expression for the shape of the buckled plate is then 



w = A^l-cos (y+By n ) 



(7) 



Substituting equation (7) and its partial derivatives into 
the energy equation (2) and reducing, yields for the criti- 
cal value of N x 

tt 2 D 



xcr ^ 3 



(8) 



where 



C = 



( L Y J_ Tl n 3 U-l) g W ! 
+ tt s L4 (an-3) 



2 {(l-n)+[(l-n)(n+l)-(n-X)jH+ (l-n)n- H s } 



TT' 



1 2 1 
— + — — H + — = — 

3 n+2 2n+l 



(9) 



Substituting for (j, in equations (6) and (9), the 
accepted value for aluminum alloys |J, = l/3 
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C = 



(L/h) 



(!) 



[" n 4 -2n s +n s 1 
L 2n - 3 J 



r 



4tt 



4 



[_3 n+ 2 



H + 



2n+l 



H" 



3tt* 



2 + (5-n) H + 



n g +n .,3"! 
3a- 1 H j 



I? 



1 + -&rr H + 



n+2 



2n+l 



H' 



(9a) 



7 " Yl 

and H = - (6 a) 

2+4n 

The expression chosen for the shape of the "buckled 
plate (equation (?)) is "based on the assumption that the 
plate will "buckle into one half wave regardless of the 
L/"b ratio. It has "been demonstrated that plates loaded 
through supported edges, and with one unloaded edge sup- 
ported and the other free (case l) , will bucklo in this 
manner (roforence 1, p. 339). There is no reason to he- 
liove that fixing the loaded edges would change this char- 
acteristic of the "buckling of such a plate. That "buckling 
will always occur in the form of one half wave means that 
the minimum value of C will occur at an infinite value 
of L/h. Referring to equation (9a), it can "be seen that 
a minimum value for 0 will occur at some "finite value of 
L/h unless the coefficient of the second term is zero. 
The only significant value of the parameter n which will 
reduce this coefficient to zero is n = 7. Equation (3a) 
then reduces to 

C = 4- + 0.406 (10) 

(L/h) 3 

Timoshenko gives an approximate expression for 0 for case 
1 (reference 1, p. 340), which involves two terms, corre- 
sponding to those in equation (l). "While this approximate" - 
expression is advocated for long plates, it also gives val- 
ues of sufficient accuracy for short plates. Tho only, sig- 
nificant difference "between equation (10) and tho corre- 
sponding equation for case 1 lies in the term involving 
(L/"b.) 3 .' For case la this term is four tiimos as large as 
for case 1. 
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Equation (s) expressed in terras of critical stress 
becomes 

•„ ' -"g si -JL- 
° cr - 12 l-v s W 

Therefore, the coefficient Z, for which values are plot- 
ted in figure 1, may be expressed 

K - 12 ° " 



Case 2a 

Loaded edges (x = 0 and x = l) fixed; one unloaded 
edge (y = 0) fixed and the other unloaded edge (y = b) 
free. 



The boundary conditions for this case are: 



[1] 


w = 


0 




when 


X 




0 


and 


when x = 1 


[2] 


3w _ 

ax " 


0 




when 


X 




0 


and 


when x = L 


[3] 


IV = 


0 




when 


7 




0 






[4] 


3 w 

ay = 


= 0 




when 


y 




0 






[5] 


3 2 w 
3y 3 


3x 3 


= 0 


when 


y 




b 






[6] 


a w 


+ (2-n) 


_ 0 

3 y3 x a 


when 


y 




b 







Assume an expression for the buckled surface of the 

form 

w = f(x) F(y) 

Since boundary conditions [53 and C 6 3 are the same as for 
case la (a free edge), equation . (3) .. is also applicable to 
this case. Boundary conditions. £3] and '[4'^] - are satisf ied 
by an expression of the form f ; " ' ' , 

w = f (x) ^1-cos 
i.e. , f(y) = (l-cos ^) (ll) 
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Substituting equation (ll) and its derivatives into equa- 
tion (3) and solving for n 

cos n = *- (12) 

2(l-n) 

for |i a 1/3 (for aluminum alloys) 

n =; 1.8343 

V7h.il e in cases 1 and la the plate can buckle into 
only one half wave, in cases- 2 and 2a the buckled config- 
uration nay contain one or more .half waves, depending on 
the proportions of the plate. In solving this problem, 
two general forms of configuration have been considered. 
First, it was assumed that the.;:plate buckled into one 
half wave. Then a solution was obtained which covers buck 
ling into any even number of half waves, the solution be- 
ing carried out specifically for two half waves. Buckling 
into a number of half waves greater than 2 was not con- 
sidered, since in such cases the effect of the condition 
of restraint at the loaded edgos is negligible, 

...Coasider ing the plate to buckle into one half wave, 
the boundary conditions [l] and [2], defining the condi- 
tions at tho loaded edges, will be satisfied by the equa- 
tion 

f(x) = A (l-cos 
and the expression for the buckled surface will be 

w = a (i_ cos ^)(i_ cos a) (13) 



Substituting equation (13) and its partial derivatives in- 
to the energy equation (2), and solving for the critical 
value for gives (for aluminum alloys) 



where 



«xcr - t> -g^- 



0 = £ + 0.0993^^ 3 + 0.624 (14) 
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As in case la 




For the type of buckling involving an even number of 
half waves « the configuration nay be approximated by the_ 
equation 

w = A (sin + 3x1 )(l- COS (15) 



vfhoro «T and B are coefficients to be evaluated from 
tho boundary conditions. In this instance the axes are 
takon as shown in figura 3, and the integrals with re- 
spect to dx in tho energy equation (2) are takon betwoon 
the Units 0 and L/2. Boundary conditions [l] and [2] 
become 

(1) w = 0 whan x = ±L/2 

(2) §1=0 whon x = ±L/3 
Pron boundary condition 1 

B = -8 sin J 

Evaluating J fron boundary condition [2], yields tho 
transcendental equation 

J - 3 tan J (16) 

The various roots of this equation give values correspond- 
ing to various even numbers of half waves of the buckled 
configuration. The lowest root corresponds to a buckled 
forn of two half waves. This root is 



J = 4.0782 



The critical value for F x can bo found as before, 
by substituting equation (15) and its partial derivatives 
into the energy equation. Such a substitution, in which 
the coefficient J has been replaced by its lowest root, 
gives for the critical load for aluminum alloys (fj, = l/3) 



where 
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0 = 5.626 + 0.0178 f-V + 0.638 (17) 
/ L\ a \ b/ 

W 

As before, the Z values shown in the curves are 
related to the G values as 



Case 4 

Loaded edge3 (x = 0 and x = L) supported; one 
unloaded edge (y = 0) fixed and the other unloaded edge 
(y = b) supported. 

The boundary conditions for this case are: 

[1] w = 0 when x = 0 and when x = L 

[2] -||f- + |o, |2x = o when x = 0 and when x = I 

[3] w = 0 when y = 0 

[4] |~ = 0 when y = 0 

[5] w = 0 when y = b ' 

r£l s g w a 2 w n 
1 Bys ^ dx 2 = 0 whea y = * 

Following the method of solution of the differential 
equation (l) described by Timoshenko (reference 1, p. 337) 
in which tho solution is taken in the form 

w = f (y) sin ^~ 

L — 

the general expression for f (y) • can be written 

f (y) = Cj. e~ ay + C a e +ay " + 0 3 cos £y + C 4 sin Py 

in which 




a / "75 — + / r-=~ a- rid B = 



m 3 TT s /-*x m 3 TT a 



is v 1,2 y D I> 2 

and G x , C 2 , C 3 , and G 4 are coefficients to be evaluated 
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from the boundary conditions. From boundary conditions 
[3] and [4] it follows that 

C 1 ~ - and Co = - 

1 2a 9 2a 



and f(y) can be written 

f(y) = A(cos£y - coshay) + B (singy - -^sinhay) 

CO 

From boundary conditions [5] and [6] two simultaneous 
equations are obtained. The critical stress can then be 
obtained by equating to zero the determinant of these equa- 
tions in A and B. This manipulation results in the 
transcendental equation 

(a 3 + £ s ) singb coshab - ^ aP + cosgb sinhab 

For given values of I»/b, this equation can be solved 
for ^ xcr » using for the number of half wavos (m) the in- 
teger which results in minimum values. The critical stress 
can be expressed in tho form 



°cr = K 



l-(j, s 



and the coefficient S thus evaluated for different ratios 
of L/b. 
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